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ABSTRACT 

A significant requirement for studies in Earth Physics is the reso- 
lution of trends from observations made over minimal time spans. 
The quantification of any surface motions with magnitudes not ex- 
ceeding a few cm per year calls for the reference of such measure- 
ments to co-ordinate systems whose behaviour is unambiguous. 
The perturbation of these systems must therefore be capable of 
accurate definition. 

A simple system of reference axes is defined for possible use in 
high precision geodetic studies over long periods of time for pro- 
grams in Earth physics. The proposed system is based on the 
gravitational and dynamic characteristics of the axis of rotation and 
the Earth's center of mass as defined instantaneously at a given 
epoch. Techniques are outlined for its continuous representation 
over time intervals of significance for studies in Earth physics. 
The relationship between the proposed system and the representa- 
tion of extra-terrestrial objects using the celestial sphere concept 
is also discussed. 
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GEODETIC REFERENCE SYSTEMS 
FOR LONG PERIOD STUDIES IN EARTH PHYSICS 


1. INTRODUCTION 

The evolution of geodetic techniques in the past two decades can be described 
without excessive exaggeration, as being spectacular. The scope of geodetic 
science whose practical applications prior to the second world war were re- 
stricted to regional studies encompassing areas no larger than a single nation, 
were extended to continental extents in the first instance. The much sought for 
connections between the various continents were subsequently achieved by the 
use of artificial Earth satellites. The precision of such determinations of posi- 
tion based on optical techniques appears to be at the level of lp.p.m. 

The co-ordinates (so established) for geodetic stations at the time of writing 
are not sufficiently sensitive for quantitative studies in Earth physics over 
periods of time shorter than 50 years. It is also uncertain whether the stations 
used are designed for unambiguous re-occupation after long periods of time. 

The application of geodetic techniques for the determination of changes in posi- 
tion for Earth physics can be summarized as follows. 

Given a system of observing stations Pj at the Earth’s surface whose relative 
positions are defined on some co-ordinate system Xj at epoch (r = t) by the co- 
ordinates Xjj on the basis of measurements made during the epoch; if the rela- 
tive disposition of the same points were re-determined as Xij at some later 
epoch (r = t + dt) and based on a co-ordinate system Xj, can the difference 

AX lj = Xij-Xjj (1) 

be interpreted in the context of relevance to the quantification of plate motion 
and crustal distortions? As estimates of gross motion (i.e. , long wave effects) 
are currently placed at 10 cm yr - ^ in sensitive areas, the quantities AXy are 
likely to be o j 1 m } if dt is approximately 1 decade. This is almost an order 
of magnitude smaller than the estimated precision of the best global geodetic 
networks available at the present time. 

It is commonly held that laser ranging techniques to appropriately equipped 
spacecraft in near Earth orbit, and to the moon should, in the first instance, 
afford ranges with a precision in excess of ±10 cm. A major geodetic task in 
the foreseeable future is the development of techniques for the efficacious use 
of the data obtained from such systems for the definition of changes in position, 
and in a manner which lends itself to use in Earth physics with the minimum 
ambiguity. 
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The definition of the station sites on ground in order that unambiguous re- 
occupation is possible at the 1-2 cm level after an extended period of time, 
poses a major problem in addition to the one outlined earlier. This question 
falls outside the scope of the present development, and it will be assumed that 
any errors on this account are less than 1 cm. A secondary problem is the 
choice of observing station sites so that changes in its co-ordinates bear a max- 
imum correlation to the motion characteristics of the crustal region it is sup- 
posed to represent. It would, for example, be undesirable if the changes AXjj 
at one of the stations in a global network but with a limited number of ground 
sites, were to be dominated by local effects with short wavelengths rather than 
those characteristics of the appropriate part of the Earth’s surface it purports 
to sample. 

The changes AXy may still reflect aspects other than those of relative motion 
even if the stability of station locations over long time spans has been clearly 
established. One major problem in four dimensional geodesy is the establish- 
ment of the relationship between the systems of reference used for each of the 
epochs of measurement. The ensuing development shows how it is possible to 
use a system which is based on certain instantaneous gravitational and dynamic 
characteristics of the Earth, as summarized in Appendix 8. 

The Earth is not a rigid body and the possibility of the constant re-distributions 
of mass should be taken into consideration, along with the associated dynamic 
consequences, when designing a system of reference for long period use. One 
of the tasks of four dimensional geodesy is to provide an unambiguous represen- 
tation of the perturbation of co-ordinate systems with time. To this end it is 
necessary to define a space for which the former provide a frame of reference. 

The system of reference described in the following sections refer to Euclidian 
space and can be related to existing concepts with minor modifications. It also 
has the advantage of being directly related to the Liouville equation (e.g. , Munk 
& MacDonald 1960, p. 10). 


2. EARTH SPACE 

All geodetic determinations to the present day have dealt with the definition of 
position on a rigid Earth. Any motion the Earth may have with respect to an 
inertial reference system in galactic space has therefore been considered to be 
effectively reversed. Earth space is defined as the Euclidian space which has 
the same rotational and galactic motion as the Earth. All geodetic co-ordinates 
established to date are Earth space co-ordinates. The effect of galactic motion 
is of relevance only when observations to extra-terrestrial objects are involved. 
The resulting consequences are listed by Mather (1972, p. 4). The various 
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aspects of galactic motion are only of indirect concern in current geodetic prac- 
tice, being absorbed as proper motion effects resulting in changes in the 
co-ordinates of the extra-terrestrial sources on the celestial sphere. 

Any geodetic determinations based on observations restricted to those between 
points on the Earth’s surface, will result in the definition of co-ordinates in 
Earth space. Observations to extra-terrestrial sources have to be corrected 
for the effect of the Earth's rotation before use for determinations in Earth 
space, unless special synchronous observing procedures have been resorted to. 
It is not difficult to show (e.g. , Ibid, p. 19) that if the Earth were a rigid body 
with every mass element rotating freely about a common axis of rotation, with 
angular velocity co, the axis should pass through the center of mass of the Earth 
(geocenter). Deviations from this concept occur due to departures of the Earth 
from a rigid body. The variation in the rate of rotation is unlikely to exceed 1 
part in 10 1 . 

Another factor which affects the definition of Earth space is the change in posi- 
tion of the axis of rotation with respect to the Earth's crust and called polar 
motion. The order of magnitude of this effect is 1 part in 10 6 . It is therefore 
convenient to model the Earth by its rigid body equivalent at any epoch in time. 
The instantaneous geocenter for this model lies on the axis of rotation. All var- 
iations are treated as perturbations which could be either periodic or secular. 
For example, the existence of departures from a rigid body model manifests 
itself as a free nutation in the motion of the pole with respect to the Earth's 
crust, with a lengthening of the period over the Eulerian equivalent (e.g. , Routh 
1905, p. 358). 

The concept of Earth space is a simple one. All the observational data available 
at the present time appears to indicate that the modeling of the Earth by its 
rigid body equivalent on lines similar to the use of a Keplerian ellipse in satel- 
lite geodesy, is an acceptable one. Deviations from this model can be treated 
as perturbations though further significant improvements in metrology may re- 
sult in a need for some multiparameter representation as being more utilitarian 
than the simple concept based on the rotation axis and geocenter. 

It could therefore be concluded that no ambiguity arises in defining Earth space 
as set out in Appendix 1 if the rotation characteristics specified are those of an 
instantaneous rigid body model of the Earth. 

Problems arise in the formulation of a suitable system of reference for the def- 
inition of position in Earth space for long period studies in Earth physics. A 
dominant concern is the extraction of meaningful information about plate motion 
and crustal distortion from equations of the type at 1 over the shortest possible 
intervals of time. The expense involved also restricts the number of locations 
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from which observations of the highest precision will be made. It would be con- 
sidered over-optimistic to put this number at even 100 by the turn of the century 
on the basis of progress being made at the present time in this direction. The 
question of a suitable system of reference must therefore be assessed in the 
context of the above factors. 

Geodetic position could be referred to an Earth space co-ordinate system de- 
fined by 

either a system of "fixed" locations at the surface of the Earth; 

or certain gravitational and dynamic characteristics of the Earth. 

The former system will require a minimum of three or possibly four stations 
held fixed as the source of the frame of reference. This is uneconomic as the 
fixed stations will thereafter make no contribution to studies in Earth physics, 
apart from defining the frame of reference. The latter system based on the 
gravitational and dynamic characteristics of the Earth will be shown to require 
only 1 station to be held fixed at the surface of the Earth to obtain a complete 
and continuous definition of position in Earth space over long periods of time. 


3. EXTRA-TERRESTRIAL REFERENCE SYSTEMS AND EARTH SPACE 

The commonly used extra-terrestrial reference systems are based on the celes- 
tial sphere concept which can be summarized as follows. The angle subtended 
at the center of the Earth by two objects whose geocentric distance tends to in- 
finity, is unaffected by either the solar orbit of the Earth or the motion of the 
objects concerned over limited periods of time. These extra-terrestrial sources 
are therefore treated as fixed in position on a spherical reference system and 
catalogued by a set of surface co-ordinates (RA, 5). The pole of the co-ordinate 
system coincides with the projection of the Earth’s pole on the celestial sphere 
at some epoch of reference (r = t 0 ). 

Observations in classical geodetic practice are oriented towards the establish- 
ment of position on a differential basis. Prior to the advent of satellite techni- 
ques, position determination was the net result of two distinct operations. The 
first involved observations as reduced to the local horizon, and called horizontal 
methods, defining as their end product, latitude and longitude, or their equiva- 
lents. The second was based on the measurement of displacements along the 
local vertical aimed at establishing the elevation of points at the surface of the 
Earth above some reference model. Horizontal survey methods give differences 
in latitude and longitude. It has been common practice to convert these differ- 
ences to "absolute" values on the basis of determinations based on the extra- 
terrestrial reference frame provided by the celestial sphere. 
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Figure 1 illustrates the basic relationships involved in transferring astronomical 
determinations into a geodetic frame of reference. If the latter were based on 
a three dimensional Cartesian co-ordinate system X i whose direction cosines 
were , assuming origin — geocenter coincidence, let P 7 be the great circle 
of zero right ascension on the celestial sphere. The diagram referred to illus- 
trates the case when the XjX 3 plane coincides with that of the great circle P 7 , 
which is equivalent to the reference meridian sidereal time being zero. 


P 



CELESTIAL EQUATOR 


(PREFERENCE MERIDIAN SIDEREAL TIME 

Figure 1. A Geodetic Reference System and the Celestial Sphere — 1. De- 
termination of Co-Latitude 0 a and Longitude X a Give 2 Angles Linking 4 
Directions; 2. 15 Direction Cosines are Linked by 12 Relations 
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Astronomical determinations of latitude (0 a ) and longitude (A a ) are obtained 
from two angles defined from observations on the celestial sphere. The first is 
the co-latitude 6 a (=-- l/2ir - 0 a ) between the pole P and the local zenith Z, GZ in 
Figure 1 having the same direction cosines £ 5i as the local vertical at the ob- 
serving station Q. On the basis of experience, it can be stated that 0 a is un- 
likely to deviate from the co-latitude 0 g on an ellipsoid which best fits the geoid 
by amounts much in excess of 0.3 mrad, the modal value being an order of mag- 
nitude smaller. This would also apply to A a which is deduced from observations 
to extra-terrestrial sources and a knowledge of the rotation characteristics of 
the Earth which are assumed to be known. The 15 direction cosines defining the 
five directions satisfy 12 equations given at 2 - 6. 


u>! 

jp 

li 

II 

k/i 

(2) 

ik ^jk = 0 (4 equations), 

(3) 

g 3i C 5i = cos0 a> 

(4) 

£ 4i £ 5i = sin0 a> 

(5) 

8 li C 4i = COS V 

(6) 


Noting that the three equations containing C 2 i can be omitted without effecting 
the development, the four conclusions summarized in Appendix 3 can be drawn. 

The relations illustrated in Figure 1 deviate from reality. Extra-terrestrial 
observations are referred to the Earth’s instantaneous axis of rotation (GPj in 
Figure 2). The quantities observed or deduced from observations are the zenith 
distance SZj and the hour angle RHA as shown in Figure 2. The co-ordinates 
(RA, 5) of the extra-terrestrial source S, on the other hand, are referred to the 
pole P r of the celestial sphere. Two problems arise on attempting to correlate 
the reference system afforded by the celestial sphere to the physical system 
characterized by the instantaneous axis of rotation. 

1. The location of the pole of the reference system in relation to the in- 
stantaneous axis of rotation . 

The problem can be viewed as the definition of the vector PjP r on the 
celestial sphere, as illustrated in Figure 3. Consequently the instan- 
taneous polar distance of S increases from (1/27T - 5) to (1/2 tt - [5 - 
A5] ). The vector PjP r can be determined in principle if the angular 
distances PjS a were known from two points S a of known position on the 
celestial sphere under the conditions normally applying to position fixes 
by intersection. 
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e= TIME INDEX ON REFERENCE MERIDIAN 

Figure 2. The Celestial Sphere, the Instantaneous Rotation System, the 
Meridian of Reference and Astronomical Determinations at 6 = 0 


2. The correlation of the "reference" meridian with the circle of zero 
right ascension . 

This can only be achieved if the concept of a meridian of reference 
were clearly defined. 

Figure 2 illustrates without unnecessary complication the problems encountered 
in positional astronomy. It has been assumed that the time index on the reference 
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Figure 3. The Change in Declination on the Instantaneous System 


meridian (0) is 0. The extra-terrestrial reference frame provides co-ordinates 
of sources with respect to a pole P r which, in the most general case, may not 
coincide with the instantaneous axis of rotation. The latter, of course, has 
measurable angular relationships on the celestial sphere to extra-terrestrial 
objects. The determination of these quantities is the observational basis of 
positional astronomy. It is conventional for P r to be the rotation pole of some 
epoch of reference (r = t 0 ). 

In the general case, P r Pj can be fully represented by two mutually perpendicular 
components 50 a , in relation to the meridian of reference (as illustrated in Fig- 
ure 2). If the problem were treated as one where it was necessary to define the 
corrections ARA to the right ascension RA and -A5 to the declination 5 of S, con- 
sider in the first instance, the question of defining the meridian of reference. 

If it were defined in terms of the axis of rotation and some point 0 at the Earth’s 
surface, it becomes obvious that the meridian of reference will keep changing 
position in Earth space as the pole moves in relation to the Earth’s crust. On 
the other hand, the plane so defined can be related to observations. 

Alternately , the reference meridian can be defined arbitrarily. The disadvant- 
age of an arbitrary definition is the difficulty in relating the reference system 
to observations. Such disadvantages could be merely conceptual rather than 
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pragmatic if the arbitrary definition were based on the Earth space location of 
the rotation axis of an earlier but recent epoch. The former alternative is il- 
' lustrated in Figure 2 which co-ordinates the great circle of zero right ascension 
to the plane of the reference meridian defined by the instantaneous axis of rota- 
tion and the zenith of O, as intersecting at O. If the co-latitude of O is 0 O with 
respect to P r , it is not difficult to show that 

ARA = - 50j cot e Q . (7) 

The correction (-AS) to the declination is given by 

AS = P i S-(l/27r-5). (8) 

Figure 3 illustrates more clearly the relation between Pi , P r and S. The eval- 
uation of AS is based on the fact that the great circle arc OS (= > P 0 ) is invariant 
between epochs. If P r were the rotation pole at an earlier epoch, as mentioned 
earlier, at which the co- latitude of the origin O was determined as 0 O , it fol- 
lows that 


PiO = 0 o (t) = 0 o (t o )-S0 2 = 0 Q - S0 2 . 


(9) 


Consideration of the spherical triangle OP r S gives 


0 o = cos -1 [cos 0 Q sin 5 + sin 0 Q cos 5 cos RA] . 


( 10 ) 


0 O is therefore a known quantity for any extra-terrestrial source represented 
on the celestial sphere, if the co-latitude of O were known at the epoch for which 
P r was the rotation pole. Further, the angle SOP r (= a), given by 


a 


sin 1 


sin RA 
sin 0 O 


cos 



(ID 


is also a known quantity. The azimuthal angle a between the instantaneous ref- 
erence meridian and the plane containing O and S changes between the epochs 
(r = t 0 ) and (r = t) by the angle 6a given by 

5a = Pj-OPj = 50j cosec (0 Q - 50 2 ) = 8(f> l cosec d Q + o {lO -10 } (12) 


if the components of polar motion between the epochs were less than 2 arcsec 
each. Consideration of spherical triangle PjOS gives 


sin(5 - A5) '= cos(0 Q - 5 0 2 ) cos 0 O + sin(0 Q - 50 2 ) sin 0 O cos(a + da). (13) 
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On expanding the above equation by Taylor’s theorem, using the assumption 
that I50o: I ^ 2 arcsec, and noting the equivalent equation from spherical tri- 

angle P r OS for 5, it follows that 

A5 cos 5 = 5</> 2 [cos0 o sin cos a- sin d Q cos 0 Q ] + 5a sin a sin d Q sin \p Q + o {(5<£ a ) 2 } 

(14) 

on use of Equation 13 and a little manipulation. 

Equations 9, 11, 12 and 15 complete the transformation of co-ordinates on the 
(RA, 5) system as defined on the celestial sphere, to equivalent quantities on 
the instantaneous system, provided the vector PjP r were known, the represen- 
tation being in terms of the orthogonal component 50 a illustrated in Figure 2. 
These components represent polar motion if P r were the rotation pole at the 
epoch of reference. 

The concepts are made more complex in practice due to 

a. the rotation of the Earth; and 

b. precession and nutation. 

The problems involved concern the relation of the celestial sphere to the physical 
system at each instant in time. They are referred to (in Appendix 5) as being 
"non-geodetic in Earth space”. It is not the intention of the present develop- 
ment to deal with such matters. Current practice is to adopt a non-complex 
model and hence linearize the system. Variations in the rate of rotation can be 
treated as parameters to be recovered in any adjustment of observations. The 
question of nutation is more complex and problems may arise in modeling the 
phenomenon to the required order of precision. 

It is significant to note that changes in the pole position can be determined with- 
out any dependence of the result on the precision of geodetic station co-ordinates. 
It has been shown that a single tracking station is capable of recovering the com- 
ponent of changes in position with respect to the instantaneous pole (Smith et al. 
1972) in the meridian of the tracking station. The analysis of such determina- 
tions at a number of stations should permit the separation of the quantities 50 a 
from other factors of significance. * 

It must therefore be concluded that the quantities 50 a can be determined without 
any dependence on the process used for the determination of geodetic position. 
Thus the XjX 3 plane defined by the instantaneous axis of rotation and the single 
station O at the surface of the Earth and its position at any epoch (r = t) can, in 
principle, be related to its location at any other epoch (r = t + dt) provided the 
sequence of changes 8<p a have been defined between epochs. An additional 
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ambiguity is introduced into the definition of the XjX 3 plane,, as summarized in 
Appendix 8, by the fact that O is defined by its zenith and not its ground location. 
Consequently, any changes in the direction of the vertical due to mass re- 
distributions will introduce an additional bias in the definition of the XjX 3 plane 
if not corrected for by gravimetric means. 


4 . A GEODETIC SYSTEM OF REFERENCE FOR LONG PERIOD STUDIES 

IN EARTH PHYSICS 

A three dimensional system of reference at any epoch (t = t) is obtained on de- 
fining an origin and a plane of reference defined by two orthogonal axes in the 
plane passing through the origin. If precession and nutation were allowed for, 
and as rigidity is implicit in the concept of rotation, the instantaneous axis of 
the Earth can be shown to pass through the geocenter. It therefore follows that 
only one point on the surface of the Earth is necessary to completely define the 
reference system. 

Appendix 8 summarizes an acceptable system of reference for geodetic measure- 
ments at any epoch (r = t), all gravitational and dynamic characteristics being 
defined by their instantaneous locations. The position of the "natural" co- 
ordinate system changes with time, despite the X jX 3 plane always passing 
through O. The locations of the co-ordinate systems providing a frame of ref- 
erence at each epoch can be inter- related in the following manner. 

1. Adopt the co-ordinate system X s defined in terms of Appendix 8 for 
some initial epoch (r = t G ) as the system of reference for long period 
studies. The origin of this Xj system coincides with the location of 
the geocenter at this same epoch. 

2 . Geodetic measurements made during some later epoch (t = t) can be 
established with respect to the reference system Xi * as described in 
Appendix 8, the X1X3 plane passing through the location of the geocenter 
for this epoch as well as the ground station O which could be a tracking 
station in the case of laser ranging systems or alternately a receiving 
antenna in the case of VLBI. 

The Xi co-ordinates so determined are related to the reference system 
values Xj by the transformation 

X[(t) ■= X gi (t) + Ojj Xj (t) ( 15 ) 
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where the elements of the transformation matrix are obtained from Figure 4 
by direct resolution as 


COS 0j COS 0 2 - 

sin 0j sin 0 2 cos 0 3 

sin 0j cos 0 2 + 
cos 0j sin 0 2 cos 0 3 

sin 0 2 sin 0 3 


- cos 0 j sin 0 2 - 
sin 0j cos 0 2 cos 0 3 

- sin 0j sin 0 2 + 
cos 0j cos 0 2 cos 0 3 


sin 0j sin 0 3 


- cos 0j sin 0 3 


cos 0 2 sin 0 3 


cos 0 3 



Figure 4. Rotating Axes in Euclidian Space, Xi =Xj - e ijkXj 00 k dt+ x g i, 
X g i = Co-ordinates of Geocenter at Epoch (t = t + dt), cj, dt = 50, , 
co 2 dt = 50 2 , u> 3 dt = 50j cot 0 O , if Celestial Sphere Maintains the Rate 

of Rotation of the Earth 


The Euler angles 0 i} illustrated in Figure 4, are a conglomeration of a series 
of increments d0j evaluated between successive epochs (r = t) and (r = t + dt). 

If Xi are the co-ordinates of a station in the geodetic network on the instantaneous 
reference system Xi summarized in Appendix 8 at epoch (r = t) and Xj those of 
the same station on the instantaneous system xj at epoch (r = t + dt), the changes 
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in position of the axes in Earth space can be completely defined by 

a. 3 rotations characterized by the angular velocities ojj about the Xi axes; 

b. the motion of the geocenter between epochs, and defined by the co- 
ordinates Xgi on the Xi system. 

Figure 4 illustrates the relations at a. above. It can be shown (Mather 1973, 
Sec. 3.2) that the rotation a>, are related to the Euler angle increments d0, by 
the first set of relations in Appendix 9. A study of Figures 2 and 4 shows that 

co a dt = 8<p a a = 1,2 (17) 

u> 3 dt = 50j cot0 Q . (18) 

In addition, it can be shown that the co-ordinates Xj of a station referred to the 
instantaneous reference system Xj at epoch (r = t + dt) are equivalent to the 
values Xj (t = t + dt) on the Xi system in Earth space, the transformation being 
defined by the equations 

Xj (t + dt) = Xj (t + dt) + ejj k cjj x k dt + Xgj (19) 

which is the standard expression for the inter-relation of rotating axes in 
Euclidian space (e.g., Mather 1972, p. 26). 

Equations 15 and 19 complete the transformation of co-ordinates between those 
referred to the instantaneous Cartesian system x 4 to values on the permanent 
reference system Xj. The differences 

AXjj(dt) = Xjj (t + dt) - Xy (t) (20) 

between the co-ordinates of the i-th station as determined at epochs (r = t) and 
(r = t + dt) would provide a measure of the change in position between epochs 
of the i-th station, which was independent of the perturbations of the co-ordinate 
systems with time. It has been assumed that the changes in position of the geo- 
center with time can be determined. 


5. THE MOTION OF THE GEOCENTER 

Both Equations 15 and 19 require a knowledge of the motion of the geocenter for 
complete definition. The geocenter is the center of mass of the solid Earth, 
oceans and the atmosphere. It is assumed that the gravitational effects of all 
extra-terrestrial sources are treated as perturbations which can be modeled 
and whose effects are removed. 
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Appendix 10 summarizes the reasons why it is considered desirable to allow for 
the motion of the geocenter. This represents a resultant motion of the Earth's 
center of mass with respect to the stations comprising the geodetic reference 
network at the surface of the Earth. Two methods are available for evaluating 
this phenomenon as summarized in Appendix 10. 

The first follows from Equation 19 when the omission of the term x g i will result 
in the existence of a first degree effect in the difference between the co-ordinates 
Xi as determined on the instantaneous system at epoch (t = t) and xj at (r = t + 
dt). Thus if Xi 'j (t + dt) were defined by 

Xj'j(t + dt) = Xjj (t + dt) + ej kg 0J k x g dt, (21) 

the recorded difference in the station co-ordinates at the i-th station is 

Axrj (dt) = Xj] (t + dt) - Xjj (t) 

= Ax ij -Xgj, (22) 

Axjj being the difference which would have been obtained had condition 2 in Ap- 
pendix 8 been satisfied. The effect of the exclusion of x g j on Axjj is harmonic 
and of first degree. This technique can only be applied to data acquisition pro- 
cedures which are directly related to the natural frame of reference specified 
in Appendix 8. Its efficacy will depend on the progress made in modeling the 
various perturbing forces affecting extra-terrestrial sources in near-Earth 
orbit. 

A second technique for evaluating x g j is from the analysis of changes in absolute 
gravity determined at the 1 jugal level, at a network of observing platforms 
spanning the globe. The method as described by Mather (1973, Sec. 3.3), is 
summarized in Appendix 11. It could be speculated that local effects may well 
be an order of magnitude larger than global effects and the success with which 
models can be defined for such local effects will be a critical factor in the 
meaningful utilization of this method. 

It may in fact be more desirable, as in astro-geodesy, to trade a slight diminu- 
tion (5 - 10 Mgal) in observing precision for a more extensive coverage from a 
less complex instrument which is more economic to operate in both time and 
costs. The required information may well be obtained by the judicious selec- 
tion of the permanent station sites in geologically inactive regions away from 
development. 

As the problems associated with each of the techniques outlined are not similar, 
it may be necessary in the first instance to use both methods simultaneously in 
order that any problems of consequence are identified and eliminated. 
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6. CONCLUSION 


The discussion in the preceding five sections has summarized the principles 
underlying the definition of the geodetic positions of a network of stations at the 
surface of the Earth without ambiguity over long periods of time by adopting the 
natural frame of reference afforded by the dynamic and gravitational character- 
istics of the Earth's rotation, and a single point at the Earth's surface. The 
sequence of development is outlined in Appendices 1 to 11. 

A question that has not been covered is the scale of the space and the possibility 
that it may change with respect to immutable standards over long periods of 
time. The implications of a variable space metric in geodetic studies has been 
discussed by Mather (1972, p. 32). Two possible effects are the following. 

1. The velocity of light c may change with respect to some immutable 
standards of length and time. 

2. Variations may occur in the gravitational constant G. 

The entire system is scaled by the instantaneous value of the velocity of light 
as both time and length standards are based on it. No inconsistency arises in 
incorporating gravity measurements into the overall context of geodetic deter- 
minations as c defines the framework for the determination of both absolute 
gravity and G. The main geodetic consequence of time variations in c is the 
occurrence of a variable set of units in defining both dflj described in Appendix 
9, and hence in the Euler angles 0 jt as well as the Earth space co-ordinates. 

In the latter case, the differences Axjj described in Equation 22 would contain 
an effect which in terms of harmonics, is of zero degree. 

The possibility of variations in G at the level of 1 part in 10 1 1 yr -1 have been 
widely discussed and a summary of its consequences on the determination of 
position is available (Ibid, p. 134). The net effect, if any, will manifest itself 
as a harmonic of zero degree in determinations of changes in absolute gravity, 
and hence x g j in Equation 19. 

It can therefore be concluded that changes in Earth space co-ordinates of geo- 
detic stations at the surface of the Earth which are of zero degree could well 
have nothing to do with changes in position, though there is the possibility of 
their occurrence as a consequence of a uniform expansion of the Earth. In view 
of uncertainty of the origin of such effects, it would be advantageous from all 
points of view to filter out zero degree changes in Earth space co-ordinates as 
a separate information matrix. The resulting Earth space co-ordinates would 
be of direct application and relevance to those studies in Earth physics concerned 
with plate motions and crustal distortions. 


15 



7. ACKNOWLEDGMENTS 


The author is grateful for the opportunity afforded by a National Academy Sci- 
ences Resident Research Associateship at the Geodynamics Program Division 
of National Aeronautics & Space Administration’s Goddard Space Flight Center, 
Greenbelt, Maryland, to be able to present this paper, while on leave of ab- 
sence from the University of New South Wales, Sydney, Australia. 


8. REFERENCES 

1. Mather, R. S., 1972. Earth Space, UNISURV Rep . G-17, pp. 1-41, Univ. 
of New South Wales, Kensington, NSW, Australia. 

2. 1973. Four Dimensional Studies in Earth Space, Bull . 

Geodes , (in press); also GSFC Preprint X-553-72-230, July 1972. 

3. Munk, W. H. and MacDonald, G. J. F., 1960. The Rotation of the Earth , 
Cambridge Univ. Press. 

4. Routh, E. J. f 1905. Advanced Dynamics of Rigid Bodies, Dover, New 
York. 

5. Smith, D. E., Kolenkiewicz, R. , Dunn, P. J., Plotkin, H. H. and Johnson, 
T. S. , 1972. Polar Motion from Laser Tracking of Artificial Satellites, 
Science, 178, 405-6. 


16 



APPENDIX 1 


EARTH SPACE 

DEFINITION 

THE EUCLIDIAN SPACE HAVING THE SAME ROTATIONAL 
AND GALACTIC MOTION AS THE EARTH 

ASSUMPTION 

EARTH IS A RIGID BODY HAVING AN AXIS OF ROTATION 

DIFFICULTIES ENCOUNTERED FOR HIGH PRECISION GEODETIC 
STUDIES 

1. DEPARTURES FROM THE RIGID BODY MODEL 

2. INSTANTANEOUS AXIS OF ROTATION HAS PERIODIC AND 
POSSIBLY SECULAR MOTIONS WITH RESPECT TO 
EARTH'S CRUST 

3. THE UNAMBIGUOUS DEFINITION OF A SYSTEM OF 
REFERENCE FOR LONG PERIOD STUDIES 
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APPENDIX 2 


GEODETIC PRACTICE TO 1 P.P.M. IS 

1. MARGINALLY SENSITIVE TO POLAR MOTION AT 0.1 ARSEC LEVEL 

2. UNAFFECTED IN A GLOBAL SENSE BY DEPARTURES OF THE 
EARTH FROM A RIGID BODY MODEL 

TIDAL EFFECTS <60 cm 
PLATE MOTIONS <10 cm yr' 1 

IF RESTRICTED TO TIME SPANS OF LESS THAN 50 YEARS 

AND IF 

a. THE ROTATION OF THE EARTH 
AND 

b. PRECESSION AND NUTATION EFFECTS ON EXTRA- 
TERRESTRIAL REFERENCE SYSTEM 

WERE ALLOWED FOR WHEN INCORPORATING INFORMATION 
DEDUCED FROM CELESTIAL SPHERE CONCEPTS INTO THE 
GEODETIC REFERENCE SYSTEM 
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APPENDIX 3 


CONCLUSIONS 

1. ASTRONOMICAL DETERMINATIONS PROVIDE 
NO SCALE 

2. THREE DIRECTIONAL PARAMETERS HAVE TO BE 
ARBITRARILY ASSIGNED TO CONVERT THE 
ASTRONOMICAL DETERMINATIONS TO 
EQUIVALENT GEODETIC QUANTITIES 

3. THE X x X 3 PLANE IS DEFINED BY THE EARTH’S 
ROTATION AXIS AND POLAR GREAT CIRCLE ON 
CELESTIAL SPHERE WITH ZERO RIGHT ASCENSION 

4. 6 *= 0 IMPLIES THAT THIS GREAT CIRCLE LIES 
WHOLLY IN THE PLANE OF THE REFERENCE 
MERIDIAN 

"PREFERENCE MERIDIAN SIDEREAL TIME 
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APPENDIX 4 


TO OBTAIN EARTH SPACE CO-ORDINATES 
FROM EXTRA-TERRESTRIAL CONSIDERATIONS 
MODELS ARE REQUIRED 

FOR 

1. EARTH ROTATION 

2. PRECESSION AND NUTATION 

IF THESE ARE ALLOWED FOR 
8 <f> a REPRESENT POLAR MOTION 

HENCE u>j ARE DEFINED 
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APPENDIX 5 

RESULTANT PROBLEMS 

NON GEODETIC IN EARTH SPACE 

1. EFFECT OF PRECESSION AND NUTATION ON THE 
VECTOR P i P r 

2. VARIATION IN THE EARTH’S RATE OF ROTATION 

GEODETIC IN EARTH SPACE 

1. POLAR MOTION 

2. MOTION OF THE OBSERVING STATIONS 

3. CHANGES IN THE DIRECTION OF LOCAL VERTICAL 
DUE TO MASS RE-DISTRIBUTIONS 

4. DEFINITION OF THE REFERENCE “MERIDIAN” 

5. MOTION OF THE GEOCENTER WITH RESPECT TO 
A GEODETIC REFERENCE NETWORK AT THE 
SURFACE OF THE EARTH 
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APPENDIX 6 


THE INSTANTANEOUS X,X 3 PLANE 

SATELLITE RANGING SYSTEMS 

1. DO NOT DIRECTLY DEPEND ON CELESTIAL SPHERE 
SYSTEM OF EXTRA-TERRESTRIAL REFERENCE 

2. XjX 3 PLANE IS UNAMBIGUOUSLY DEFINED IF 0 IS A 
STATION IN THE GEODETIC NETWORK AND VARIATIONS 
IN THE RATE OF ROTATION ARE ALLOWED FOR WHEN 
ADJUSTING OBSERVATIONS 

3. ALL OBSERVATIONS REFERRED TO INSTANTANEOUS 
GEOCENTER 

ASTRONOMICAL SYSTEMS 

AN ADDITIONAL AMBIGUITY IS INTRODUCED BY 
CHANGES IN THE DIRECTION OF THE VERTICAL AT 0 
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APPENDIX 7 

OPTIMUM CHARACTERISTICS 
IN A GEODETIC REFERENCE SYSTEM 
FOR LONG PERIOD STUDIES IN EARTH PHYSICS 

1. OBSERVATIONS CAN BE RELATED DIRECTLY TO THE 
REFERENCE SYSTEM 

2. A MINIMUM NUMBER OF POINTS SHOULD BE HELD FIXED 
AT THE SURFACE OF THE EARTH 

- PREFERABLY ONE 

3. STATIONS DEFINING THE REFERENCE SYSTEM SHOULD 
HAVE A PERMANENT CHARACTER 

- GEODETIC OBSERVATORIES 

4. THE REFERENCE SYSTEM ADOPTED SHOULD BE SUCH THAT 
CHANGES IN CO-ORDINATES BETWEEN EPOCHS CAN BE 
SIMPLY RELATED TO SECULAR VARIATIONS IN POSITION 
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APPENDIX 8 


AN ACCEPTABLE SYSTEM IS OBTAINED AS FOLLOWS 

1. AT ANY EPOCH (r=t) REPLACE EARTH BY AN 
INSTANTANEOUS RIGID BODY MODEL 

2. SYSTEM OF REFERENCE IS THREE-DIMENSIONAL 
CARTESIAN CENTERED AT GEOCENTER 

3. X 3 AXIS COINCIDES WITH AXIS OF ROTATION 

4. XjX 3 PLANE IS DEFINED BY THE AXIS OF ROTATION 
AND ONE LOCATION AT THE EARTH’S SURFACE 
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APPENDIX 9 


GOAL OF GEODETIC REFERENCE SYSTEM 
TO PROVIDE SIMPLE INTER-RELATION OF CO-ORDINATES 
BETWEEN WIDELY SEPARATED EPOCHS (r=t 0 ) AND (r=t) 


ACHIEVE BY 

1. USING A SYSTEM OF EULER ANGLES 0j 

2. DEFINING THE MOTION OF THE GEOCENTER 


EULER 

ANGLE 

INCREMENTS 


ft T = o>| CO 2 OJ 3 | 

d 0 T = d 0, d 9 z dt 
sin 0 2 sin 
R = cos 0 2 sin 
cos 9j. 


0 cos 0 2 

0 -sin 0 2 

1 0 


CONVERSION 

OF 

INSTANTANEOUS 

CO-ORDINATES 

TO 

REFERENCE 

Xi 

SYSTEM 



Xj (t) = Xgj(t) + aijXj(t) 

ct j j ARE ELEMENTS OF ARRAY 

cos 0 | cos 9 2 - - cos 9 ( sin 9 Z - 

sin 9 1 sin 9 2 cos 9 3 sin 9 ( cos 0 2 cos 

sin cos 0 2 + -sin si n + 

cos 9\ sin 02 cos 0 3 cos 0| cos 02 cos 03 

sin 02 sin ©3 cos 0 2 sin 0 3 


sin 0 | sin 0 3 


- cos 0 | sin 0 3 


cos 0 3 
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APPENDIX 10 


THE MOTION OF THE GEOCENTER 

IF MOTION IS IGNORED 

1. CO-ORDINATE CHANGES OF SURFACE GEODETIC NETWORKS 
CONTAIN INFORMATION ABOUT MOTION OF GEOCENTER, 

IF OF SIGNIFICANT MAGNITUDE 

2. THIS INFORMATION IS BETTER EXCLUDED IF CO-ORDS ARE 
USED FOR REGIONAL EARTH PHYSICS STUDIES 

TECHNIQUES FOR EVALUATION IN THE FUTURE 

1. ANALYSIS OF CHANGES IN STATION CO-ORDINATES BETWEEN 
EPOCHS FOR FIRST DEGREE HARMONIC EFFECTS 

2. ANALYSIS OF CHANGES IN ABSOLUTE GRAVITY AT THE 
±1/*GAL LEVEL BETWEEN EPOCHS FOR FIRST DEGREE 
HARMONIC EFFECTS 
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APPENDIX II 


AT 

EPOCH r = 

t 




AT i-th STATION 




(Agj) t 

A 

= Qi + 

A 

c i 

A 

Xi 

4 


1 

GRAVITY 

1 

OBSERVED 

1 

ELEVATION 

T 

NORMAL 


ANOMALY 

GRAVITY 

EFFECT 

GRAVITY 

AT 

EPOCH r = 

t+dt 




( Ag 

- (Qj + Sgj) + 

Cj - Xi 

i i 



t 

GRAVITY 

ANOMALY 

t 

OBSERVED 

GRAVITY 

t 

AS ABOVE 




= ( Ag j) t + Sg j 




Sgi DUE TO 

1. MOTION WITH RESPECT TO ROTATION AXIS 

1/xgal ~ lm 

2. RADIAL MOTION RELATIVE TO GEOCENTER 

1/xgal ~ 3mm 

3. LOCAL EFFECTS e.g. INCREASE IN WATER TABLE 

1/u.gal ~ 2-3cm 

GLOBAL ANALYSIS FOR A FIRST DEGREE HARMONIC AT 
(2) GIVES MOTION OF GEOCENTER BETWEEN EPOCHS 
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APPENDIX 12 


SCALE 

ALL MODERN GEODETIC SYSTEMS ARE 
SCALED BY THE VELOCITY OF LIGHT C 

CONCLUSION 

IT IS IMPORTANT THAT ZERO DEGREE 
CHANGES IN CO-ORDINATES BE 
REMOVED 
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APPENDIX 13 


ADVANTAGES OF SYSTEM PROPOSED 

1. UNAMBIGUOUS IF SCALE IS SPECIFIED 

2. NOT COMPLEX 

3. DIRECTLY RELATED TO OBSERVATIONS 

4. NO SIGNIFICANT CHANGES IN CONCEPT FROM 
SYSTEMS CURRENTLY IN USE 

5. DIRECTLY RELATED TO LIOUVILLE EQUATION 


NASA— GSFC 
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